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Abstract 

This paper presents a distributed algorithm to simultaneously compute the di- 
ameter, radius and node eccentricity in all nodes of a synchronous network. Such 
topological information may be useful as input to configure other algorithms. Pre- 
vious approaches have been modular, progressing in sequential phases using building 
blocks such as BFS tree construction, thus incurring longer executions than strictly 
required. We present an algorithm that, by timely propagation of available estima- 
tions, achieves a faster convergence to the correct values. We show local criteria for 
detecting convergence in each node. The algorithm avoids the creation of BFS trees 
(~*\ and simply manipulates sets of node ids and hop counts. For the worst scenario of 

^ variable start times, each node i with eccentricity ecc(?) can compute: the node eccen- 

O tricity in diam(G) + ecc(i) + 2 rounds; the diameter in 2 diam(G) + ecc(i) + 2 rounds; 

and the radius in diam(G) + ecc(i) + 2 radius(G) rounds. 
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1 Introduction 

oo 

This paper presents a distributed algorithm to simultaneously compute the diameter D, 
radius R and node eccentricity ecc(i) in all nodes of a network. An early knowledge of this 
topological information is useful since it is often used as input to other algorithms. For 

t-h instance, the diameter or eccentricity can be used to simplify termination in leader election 

algorithms [8j and calibrate time-to-live parameters [7]; the radius and eccentricity allow 
determining center nodes [5], which are nice candidates to serve as coordinators in other 

^ distributed algorithms. 

We assume a synchronous network model, while allowing variable start times, in which 
one or more nodes can start the algorithm with no prior coordination. The algorithm 
is designed to be fast in a precise sense; we are concerned with, not just asymptotic 
complexity, but exact bounds in the number of rounds. 

The classic approach to this problem |8j is to compute the eccentricities by parallel 
construction of breadth first search (BFS) trees rooted at each node. Once eccentricities 
are known, each BFS tree can be reused to do a global computation, starting from the leafs 
and converging to each root node, allowing each to compute the maximum and minimum 
eccentricity (the network diameter D and radius R). Considering a graph G = (V,E), 
this classic approach has total message complexity of 0(|V| \E\ log |V|) bits. The diameter 
and radius are known at all nodes in at most 4D + 2 rounds. 
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These time bounds can be improved if one departs from this modular multi-phase 
approach, where BFS trees are first constructed to compute eccentricities. This paper 
introduces an algorithm that propagates candidate values in a timely and continuous 
fashion, resulting in a faster convergence to the correct values. The challenge in this 
strategy is that a suitable termination method must be devised to detect, in each node, 
when the candidate values have converged. Under the same message complexity of the 
classic approach, the proposed algorithm reduces the number of rounds to compute the 
diameter to at most 3D + 1 rounds and the radius to at most 2D + 2R rounds. To be 
more precise, with this algorithm each node i with eccentricity ecc(i) computes: 

• the node eccentricity at most by round diam(G) + ecc(i) + 2; 

• the diameter at most by round 2 diam(G) + ecc(i) + 2; and 

• the radius at most by round diam(G) + ecc(z) + 2radius(G). 

The paper is organized as follows. Section [2] presents the computing model and intro- 
duces notation. The algorithm is presented in Section [3j Example runs of the algorithm, 
proofs of local convergence criteria and global convergence bounds are also included in 
this section. The related work is discussed in Section |4j and the conclusions are presented 
in Section [5j 

2 Network Model and Notation 

We assume a synchronous network model, similar to the one described in [8]. The network 
is composed by a set of nodes connected by links, which we assume to be bidirectional; i.e., 
we have a simple, connected, unweighted, undirected graph G = (V, E) with \ V\ > 2 nodes 
and \E\ > 1 links. We assume globally unique identifiers for nodes, but no knowledge of 
the network topology or the number of nodes. 

Computation proceeds in synchronous rounds. At each round, nodes first look at their 
state and compute what messages are sent, through a message- generation function; then 
nodes look at their state and messages received and compute the new state after the round, 
through a state-transition function. We assume no link or process failures. In order to 
obtain an asynchronous version of the algorithm a synchronizer a p] can be used. We 
operate under a maximum bandwidth of 0(|U| log \V\) bits, per link per round. 

We assume the general case with variable start times. Nodes start as quiescent, a state 
in which they do not send messages nor transition to different states. Nodes wakeup when 
they receive a message from a special environment node (not part of G, connected to every 
node), or from an already active node. 

We use d(i,j) to denote the distance between nodes i and j (the length of the shortest 
path between nodes i and j); ecc(i) for the eccentricity of node i (the maximum d(i,j) 
between node i and any other node j); diam(G) for the network diameter (the maximum 
eccentricity over all nodes); radius(G) for the network radius (the minimum eccentricity 
over all nodes); and nbrs(i) for the set of neighbors of node i (nodes connected to node i 
by a link). 
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3 Algorithm 



The algorithm is presented in Figure [TJ At each round, a node i sends the same message 
to all its neighbors (state variable Oj). A message is a non-empty set of tuples; the empty 
set represents absence of a message. The tuples in a message can be (bfs, _, _), (diam, _) 
or (rad,_), where bfs, diam and RAD are constants. Nodes do not need to distinguish 
between messages that arrive from different neighbors; the second parameter of the state- 
transition function (parameter Mi) is the set of messages received by node i from all 
neighbors. 

Each node when awoken broadcasts a BFS message with its id and a hop counter, which 
starts at 0. Nodes keep the set of ids of all received BFS messages. When a node receives a 
bfs message from a node not yet known, it increments the hop counter and rebroadcasts 
it. 

Nodes know their eccentricity is at least the largest hop count received in a bfs mes- 
sage, which they keep in a variable (ej). Nodes also keep in two other variables (di, rj) 
a lower bound estimate for the diameter and an upper bound estimate for the radius. 
When a node increases the diameter estimate; it broadcasts a diam message with the new 
value. A node increases its diameter value when (1) its eccentricity surpasses its diameter 
estimate or (2) it receives a DIAM message whose value is higher that its diameter value. 
Estimation of the radius is also driven by eccentricity, but must be deferred until each 
node detects its correct eccentricity. When that happens, and also when a lower estimate 
for the radius is received, a RAD message is broadcast. 

It is easy to see that at some point every node will have awakened; later everyone 
will have received a bfs from everyone else and will have their eccentricity stored in the 
respective variable; and later still nodes will receive some DIAM message with the network 
diameter, originating from some maximum eccentricity node (a periphery node). Similarly, 
a RAD message originating from a minimum eccentricity node (a center node) will arrive 
eventually at all nodes. 

The relevant question is convergence detection, i.e., when will nodes know that their 
eccentricity, diameter and radius variables have converged to the correct values. For this 
purpose, and inspired by the approach in [9j, nodes have a variable which stores the 
number of consecutive rounds for which no new BFS messages arrived. Later, we will show 
how this variable can be used for convergence detection. 

In order to analyze the communication complexity of the whole execution of the algo- 
rithm, we can first observe that each bfs message can be encoded in ©(log |V|) bits, since 
its dominated by the size of the ids. Each node retransmits exactly one bfs message for 
every other node, totaling 0(|V| \E\ log | V|) bits. Since the diameter can only increase at 
most D times, each node broadcasts at most D diam messages, totaling 0{D \E\ \ogD) 
bits. Similarly for RAD messages. Thus, total message complexity is 0(|V| |£'|log|V|) 
bits. 

3.1 Example Runs 

Prior to the formal proofs of the algorithm properties, we now convey some intuition by 
illustrating its execution in two different graphs. The first graph is a path with eleven 
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state variables: 

ej, node eccentricity, initially ei = 
di, network diameter, initially dj = 
rj, network radius, initially r; L = oo 
Si, status, initially Sj = quiescent 
/j, set of node ids, initially h = {} 

Cj, consecutive rounds with no new BFS, initially Cj = 
Oi, message to be sent, initially Oi = {} 

message-generation function: 

msg i ((ei,di,r i ,Si,Ii,Ci,Oi},j) = Oi j G nbrs(i) 

state-transition function: 

transi((e i ,di,ri,Si,/i,c i ,Oi),M i ) = (ej, d$, rj, s-, I[, cj, O-) 

where 

M = U{™ I m G MJ 

if S j = quiescent A M = {} then 

! ^ii r ii s ii -^it c it^i) = ( e 'h di, Ti, Si, Ii, Ci, Oi) 

else 

M' = {(BFS,i,/l + l) | (BFS,J,/l) G M,jgli} 
M" = M' U {(BFS, i, 0) | Si = quiescent} 
if M" = {} then 

I C^ = Q + 1 

else 

L^ = o 

e- = max({ei} U {h | (BFS, _, /i) G M'}) 

4 = max({di} U {d | (DIAM, d) G M} U {e'J) 

r- = min({ri} U {r | (rad, r) G M} U {e- | c- = 2}) 

s'j = ACTIVE 

I' i= IiU{j\ <BFS,j,_>GM"} 
M d = {(DIAM,^) j dj > di} 

M r = {(rad, r^) | < n} 
_ 0[ = M" \JM d UM r 



Figure 1: Algorithm. 
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Figure 2: Path graph. 




nodes, depicted in Figure [2} Nodes and 10 and are the only two nodes in the periphery, 
thus defining the diameter; node 5 is the single node in the graph center. We consider a 
run where node is activated and we will observe how the local variables evolve in nodes 
0, 5 and 10 (Figures [3j [4] and [5j respectively). 

When probing the variables at node we can observe that each two rounds a new 
BFS is received until BFS messages from all nodes have arrived. Each time a new BFS 
is received the Cj counter is reset. Thus, when a reaches two the node knows that its 
eccentricity variable ej has reached the correct final value. This happens in nodes 0, 5 and 
10 at rounds 22, 17 and 22, respectively. 

Once eccentricity is stabilized nodes start disseminating rad messages in order to 
detect the minimum eccentricity. Since local radius estimates can only decrease, one must 
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Figure 6: T shaped graph. 

establish a termination condition so that each node knows that it reached the correct 
radius. We will show that this is safely achieved when Cj > 2rj. The run in Figure [5] shows 
that no sooner than this has the radius reached its final value of 5 in node 10. 

Diameter dissemination, via diam messages, starts even if nodes are still updating 
their (monotonically increasing) estimates for eccentricity. For instance, in Figure [4] we 
see that node 5 is activated at round 5 (receiving bfs messages from nodes to 4) and sets 
both a and d% to 5 (its received hop distance from node 0, the furthest away). As nodes 6 
to 10 are activated, in the next rounds, the diameter estimate di at node 5 gets updated 
each two rounds. We will show that in order for a node to know that the diameter has 
reached its final value it needs to locally observe that Cj > 2 and Cj > d{. This happens in 
nodes 0, 5 and 10 at rounds 31, 26 and 21, respectively. One can notice that in this path 
graph the diameter is stable even before those rounds. 

However, it is easy to construct a T shaped graph where detection cannot occur before 
the mentioned round. In Figure [6] we show such a graph. Here, a critical case happens 
when a node not in a diameter defining path (for example, node 14) is the first to be 
activated. As seen in Figure [7] this leads to a run where the diameter estimation seems 
to be stable for a large number of rounds (between round 18 and 27), finally increases at 
round 28, with convergence being detected only when a > di at round 29. In the next 
section we prove that these convergence criteria are correct for general graphs. 
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3.2 Local Convergence Criteria 

We now establish results that allow a node to know, using local information, that the vari- 
ables estimating node eccentricity, network diameter and network radius have converged 
to the correct values. In the following we will use i, j, u, and v to range over node ids and 
r, n and k to range over non-negative integers. 

Definition 3.1. Aj(r) = {j \ d(i,j) < r A Sj(r - d(i,j)) = ACTIVE}. 

Ai(r) denotes the "area of visibility" of node i afteiQ round r: the set of nodes whose 
bfs messages arrive at i no later than round r. It is easy to see that Ai is monotonic: 
Ai(r)QAi(r + l). 

Lemma 3.1. After any round r , h{r) = Ai{r). 

Proof. If j G Ii(r) then i must have received no later than round r a BFS message starting 
at j; it traveled d(i,j) hops, which means that d(i, j) < r and j was active after round 
r — d(i,j); therefore j £ Ai(r). If j £ A{(r), then d(i, j) < r and j was active after round 
r — d(i,j). This implies that j's bfs message arrives at i not later than round r; therefore, 
j G h{r). □ 

Lemma 3.2. If Ci(r + 2) > 2, then Ai(r) = V. 

Proof. If Ai{r) ^ V then there are two nodes, u 6 Ai(r) and v ^4j(r) which are adjacent, 
i.e., d(u, v) = 1. This means that s u (r— d(i, u)) = active and s v (r— d(i, v)) = quiescent. 
Since u and v are adjacent, then s v (r+l— d(i, u)) = active. There are three possible cases: 
(1) d(i, u) = d(i, v), in which case v became active in round r+1 — d(i, v), and i receives the 
BFS from v at round r+1; (2) d(i, v) = d(i, u) + l, in which case s v {r+2 — d(i, v)) = ACTIVE, 
and i receives the BFS from v at either round r + 1 or round r + 2; (3) d(i, v) = d(i, u) — 1, 
cannot happen, as it contradicts s v (r — d(i,v)) = QUIESCENT. In any case, Cj(r + 1) = 

1 When referring to state variables, we do not use the expression "at round r" to avoid ambiguity 
between beginning or end of round. Throughout the paper we use "after round r" as a shorthand for 
"when round r has finished", i.e., "at the end of round r". 
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or Cj(r + 2) = 0. Therefore, since q can only increase 1 unit per round, it follows that if 
Ci{r + 2) > 2, then Ai(r) = V. □ 

Lemma 3.3. After any round r, ej(r) = max({0} U {d(i, j) \ j E Ij(r)}). 

Proof. Trivial induction on the number of rounds. □ 

Theorem 3.1 (eccentricity convergence). If Cj(r + 2) > 2, then ej(r) = ecc(i). 

Proof. Combine the previous three lemmas. □ 

Lemma 3.4. If for some r and n > 2, Ci(r + n) = n, then for all k, Ci(r + k) = k. 

Proof. After the first round r such that Ci(r + 2) = 2, by the first two lemmas Ii(r) = V . 
This implies that for r 1 > r, Ci(r' + 1) = Cj(r') + 1. □ 

Theorem 3.2 (diameter convergence). When Cj(r) > 2 and Cj(r) > di(r), then di(r) = 
diam(G). 

Proof. By contradiction. Assume c»(r) > 2 and Cj(r) > di(r) but di(r) < diam(G). By 



Theorem 3.1 ej(r) = ecc(i). Also, it is trivial that di(r) > ej(r). Then, as in a graph 
all the eccentricities between the radius and the diameter are present [3], there exists 
two nodes u and v with d(u,v) = d{(r) + 1. Assume without loss of generality that 
d(i,u) > d(i,v). From the previous lemma, for r' = r — Cj(r), it follows that Cj(r' + k) = k. 



As Ci(r' + 2) = 2, by Lemma 3.2, Ai(r') = V; therefore s u (r' — d(i,u)) = ACTIVE. 
Then, d v (r' — d(i,u) + d(u, v)) > d(u,v) (bfs from u has reached v). Furthermore, 
di(r' — d(i,u) + d(u, v) + d(i,v)) > d(u,v) (diam message from v has reached i). Since 
d(i,u) > d(i,v), then di(r' + d(n, v)) > d(u,v). Recall that d(u, v) = d, L {r) + 1, and let 
r" = r' + d(u, v) = r — c«(r) + di(r) + 1. From the assumption Cj(r) > d«(r), it means that 
r" < r, which together with di(r") > di(r) + 1 contradicts the monotonicity of d%. □ 

Theorem 3.3 (radius convergence). When Cj(r) > 2rj(r), i/ien rj(r) = radius(G). 

Proof. We assume networks with at least one link and two nodes, which means rj(r) > 1. 



If Ci(r) > 2rj(r), we have Cj(r) > 2, which means that, from Lemma 3.4, all bfss have 
already reached node i after round r' = r — Ci(r), and from r' on we have Q(r' + k) = k. 
Assume, by contradiction, that r,(r) > radius(G). Then, at most after round r'+rj(r) — 1, 
all bfss have reached some node u in the center of the network. At most two rounds later, 
after round r" = r' + ri(r) + 1, we have c u {r") > 2 and u has sent a rad message with the 
network radius. At most rj(r) — 1 rounds later, after round r'" = r' + 2r^(r) this message 
arrives at i and rj(r w ) = radius(G). But assuming c,(r) > 2r^(r) it means that r'" < r, 
which contradicts rj being monotonically decreasing. As rj results from some eccentricity 
and is always an upper bound of the radius, we must have rj(r) = radius(G). □ 
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3.3 Convergence Bounds 

We now determine upper bounds on the number of rounds for convergence of eccentricity, 
diameter and radius. Given that we have described the algorithm for the general case of 
variable starting times, what matters is the number of rounds after the first activation; i.e., 
ignoring an initial sequence of rounds with all nodes inactive. Therefore, in this section 
we consider that the first node became active after round 0; round 1 is when the first 
non-environment message is sent. 

Proposition 3.1 (eccentricity bound). Node i can determine its eccentricity at most in 
diam(G) + ecc(i) + 2 rounds. 

Proof. After round diam(G) all nodes are active, so the last BFS arrives at i at most after 



round diam(G) + ecc(i). Two rounds later the q variable reaches 2 and from Lemma 3.1 



the eccentricity has already converged. □ 

Proposition 3.2 (diameter bound). Node i can determine the network diameter at most 
in 2diam(G) + ecc(i) + 1 rounds. 

Proof. After round diam(G) all nodes are active, so the last bfs arrives at i at most after 
round diam(G) + ecc(i). Subsequently q starts increasing and after further diam(G) + 1 
rounds the local condition q > di is met. As we are considering networks with at least 



one link, i.e., diam(G) > 1, then at this round we have also > 2 and from Theorem 3.2 



the diameter has converged. □ 

Proposition 3.3 (radius bound). Node i can determine the network radius at most in 
diam(G) + ecc(i) + 2radius(G) rounds. 

Proof. After round diam(G) all nodes are active, so the last BFS arrives at i at most after 
round diam(G) + ecc(z); afterwards Cj starts increasing and after round r = diam(G) + 
ecc(i) + 2radius(G) we have Cj(r) > 2radius(G). Also, after at most round diam(G) + 
radius(G) all BFS have arrived at all center nodes; two rounds later, at most after round 
r' = diam(G) +radius(G) + 2, each center node sends a RAD message containing radius(G). 
There are two possibilities: (1) ecc(i) > 1, the RAD message from a center node j arrives 
at i at most radius(G) rounds later, which means that at most after round r" = diam(G) + 
2radius(G) + 2 we have rj(r") = radius(G); given that r" < r, then r-j(r) = radius(G), 
and the local radius convergence criteria Cj(r) > 2r^(r) is met; (2) ecc(i) = 1, in which 
case radius(G) = 1, i is a center node, and at round r' we have r,i{r') = radius(G); as in 
this case r' = r, we have Cj(r) > 2rj(r) as well. □ 

Corollary 3.1. All nodes know: their eccentricity at most in 2D + 2 rounds; the diameter 
at most in 3D + 1 rounds; and the radius at most in 2D + 2R rounds. 

3.4 Termination 

To keep the presentation clear and avoid cluttering, we did not include in the algorithm 
the mechanics of termination; i.e., each node reaching a "terminated" state in which it 
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stops sending messages. In general distributed termination is independent from reaching 
some result, and nodes may have to keep propagating messages for some time. 

In this case, however, it is easy to see that when a node has determined both the 
radius and diameter through the local convergence criteria, all neighbors will have the 
same criteria met after at most one more round. (After one more round, each node j 
neighbor from i, will have Cj with at least the same value node Cj had, and both rj and 
dj will have the same values as in node i.) Therefore, after having met both criteria for 
radius and diameter, a node needs only execute one more round and stop. 

3.5 Improving Storage Requirements 

In the previously described algorithm each node accumulates in the / variable all ids re- 
ceived in all previous rounds. Although it has made the description intuitive and stream- 
lined proofs, it means that, regardless of network topology, by the end of the execution 
each node will need to store 0(|V|) ids. 

Here we show that it is enough to keep in the state only the ids received in the two 
previous rounds. While this modification does not change the worst case space require- 
ment complexity (it still remains 0(|V|) ids for general graphs and uncoordinated start 
times), it may be useful in practice. As an example, for 2D geometrical networks (e.g. a 
geographically spread sensor network with links according to inter- node distance), under 
synchronized start times, the number of ids that arrive in a single round (and need to be 
stored) will be V|). Notice that these specific configurations also reduce the required 
channel bandwidth, and that in other specific graph topologies, these uppers bounds on 
stored state can be even more tight. 

The modification to the algorithm is trivial and consists of replacing state variable / 
by a pair I, J used as a sliding window; replacing the test j ij with j Ii U Jj and 
replacing I \ = Jj U {j | (bfs, j, _) G M"} with I\ = J; and J[ = {j | (bfs, j, _) G M"}. The 
modification is possible due to the following property of the original algorithm. 

Proposition 3.4. A node j can only receive BFS messages (bfs,«, _), originated in a 
node i activated at round r, in rounds r + d(i, j), r + d(i,j) + 1 and r + d(i,j) + 2. 

Proof. The first round, where node j can receive a (bfs, i, _) message, is r' = r + d(i, j), by 
the shortest path from i; j rebroadcasts it and at round r' + 1 it arrives at all neighbors, 
if any. Also at round r' + 1 node j may receive such a message if there exists a neighbor 
node u at the same distance from i (i.e., d(i,u) = d(i,j)); this neighbor, similarly to 
j, has received such a message at round r' and has rebroadcasted it. At round r' + 2 
node j can receive such a message if it has a neighbor u one hop further away from i 
(i.e., d(i,u) = d(i,j) + 1); this neighbor has received the message at round r' + 1 and has 
rebroadcasted it at round r' + 2. Because any neighbor u of j has stored i in the respective 
variable after at most round r' + 1, it will not rebroadcast any (bfs,z,_) message in 
any round later than r' + 2, from which such messages cannot reach j later than round 
r' + 2. □ 
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4 Related Work 



As mentioned in the introduction, our algorithm improves the classic modular approach [8], 
where BFS trees are first computed at each node and later reused to perform a global com- 



putation of the network radius and diameter (in at most AD + 2 rounds). By Corollary 3.1 



we achieve a speedup of D rounds for computing the diameter, and since D < 2R < 2D 
the speedup for computing the radius varies between 2 and D + 2 rounds. The maximum 
speedup occurs, for example, in path graphs. This improvement is achieved with the same 
message complexity of 0(|F| |.E|log|V|) bits, and the same space complexity of 0(|V|) 
ids and computation complexity per round of 0(|^| 2 ) at each node. 

The work in |10j computes the diameter under the more restrictive synchronized start 
time model, where all nodes are activated at the first round. This is a fast algorithm since 
they also disseminate candidate values for eccentricities before they converge. However, 
since they assume that all nodes are active in the first round the local termination criteria 
is much simpler. When restricted to a setting where all nodes are active in the first round 
our algorithm outputs the diameter in the same time bound. Even though we are more 
general, we have significant improvements in message and space complexity. 

Related to the computation of the radius and the eccentricity is finding the network 
center. To the best of our knowledge, within similar bounds for space and processing com- 
plexity per round, the fastest algorithm so far to find network center nodes was proposed 
by Korach, Rotem and Santoro [6]. This algorithm builds on a simpler algorithm to find a 
center of a tree, and the observation that the center of a general network must also be the 
center of its own BFS tree: the initiator node triggers BFS generation at all nodes, picks 
a center among all candidates which are centers of their own BFS, and later disseminates 
this information to all nodes. The closest center node c to the initiator node i will receive 
the confirmation that it is indeed a center of the network by round 4 radius(G) + d(i, c) + 1, 
with total message complexity of 0(|V| \E\ log | V|) bits. Notice that, in contrast to this, 
our algorithm is fully symmetrical, not requiring a distinguished node as initiator (which 
would have to be chosen in some way, e.g. by a leader election). Even discounting this 
factor, our more general algorithm has no time penalty: in fact, it can be shown that it 
even improves this upper bound by at least one round. 

A more general problem than computing the eccentricities, radius and diameter is 
that of computing the distance matrix or all-pairs shortest path matrix in a network. In 
fact, relying on this connection, a faster distributed variant of our algorithm could be 
designed to compute the radius and diameter: instead of propagating just the distances, 
each node propagates sets of neighbors - at most by round 2 diam(G) + 2 all nodes would 
be able to determine the full topology of the network and run a standard unweighted all- 
pairs shortest path algorithm with computation complexity 0(|V| |-E|) at the last round. 
Unfortunately, even assuming that computation complexity per round is negligible, this 
algorithm requires 0(|y| 2 ) space complexity at each node, which is impractical for large 
networks. Moreover the message complexity increases to e(\E\ z log\V\) bits, and requires 
a bandwidth of 0(|-E| log | V|) bits per link. 

With such large bandwidth it is possible to decrease the message complexity using 
more elaborated approaches. For example, Kanchi and Vineyard [5] propose a distributed 
algorithm to compute the all-pairs shortest path matrix with message complexity of 
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0(|V| |£'|log|V|) bits: a spanning tree is first computed using the algorithm proposed 
by Awerbuch |2j and later reused to propagate the topological information to a root node 
that computes the distance matrix and disseminates it to all nodes. The tradeoff for this 
optimization is, unfortunately, a substantial increase in the number of rounds, although 
still 0(\V\). 

Message complexity can also be reduced by trading off accuracy. For example, Gu 
and Cheng [4| propose a distributed algorithm to compute an estimate diam(G) for the 
diameter, such that diam(G) < diam(G) < diam(G) + 2. Unfortunately, the bandwidth 
requirements are slightly worse than ours and the improvement in message complexity 
does not extend to time complexity: although the authors do not quantify this measure, 
it is clear from the algorithm presentation that the number of rounds until completion 
is substantially larger than ours. Moreover, this algorithm also requires a distinguished 
node as initiator. 

5 Conclusions 

In this paper we propose a time efficient algorithm that computes in all nodes of a network 
the values of the node eccentricity, and the network diameter and radius. The algorithm is 
very flexible in the sense that it does not require a distinguished node, one or more nodes 
can initiate the computation, and concurrent initiations have no detrimental impacts on 
the various bounds and complexities. 

Under the same communication, space, and computation complexity, the presented 
algorithm significantly improves existing time bounds for the diameter and radius com- 
putation. It also slightly improves the time bounds for the special case of finding center 
nodes, while relaxing the need for a special initiator. 

The key to the improvement was to abandon the traditional modular approach. In- 
stead, our algorithm relies on a very early propagation and aggregation of the maximum 
and minimum candidate eccentricities, even before these values have stabilized. Together 
with adequate convergence detection criteria, this allowed a simple and fast approach to 
the computation of these distances. 
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